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The Fortuin-Kasteleyn and heat-bath damage spreading temperatures Tfk{p) and Tdsip) are 
studied on random bond Ising models of dimension two to five and as functions of the ferromagnetic 
interaction probability p; the conjecture that Tds{p) ~ Tfk{p) is tested. It follows from a statement 
by Nishimori that in any such system exact coordinates can be given for the intersection point 
between the Fortuin-Kasteleyn Tfk{p) transition line and the Nishimori line, [pnl,fk ,Tnl.fk]- 
There are no finite size corrections for this intersection point. In dimension three, at the intersection 
concentration [pnl.fk] the damage spreading Tds(p) is found to be equal to Tfk{p) to within 0.1%. 
For the other dimensions however Tdsip) is observed to be systematically a few percent lower than 
Tfk{p). 
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INTRODUCTION 



The physics of Ising spin glasses (ISGs) in the "or- 
dered" regime below the freezing temperature Tg has 
been intensively investigated for decades; the param- 
agnetic regime above Tg has attracted less attention. 
However, in addition to the standard ordering (or freez- 
ing) transitions other significant "critical" temperatures 
within the paramagnetic regime can be defined opera- 
tionally and estimated numerically with high precision. 

We have studied random-bond Ising models (RBIM) in 
dimensions two, three, four and five having ferromagnetic 
near neighbor interactions with probability p and antifer- 
romagnetic interactions with probability 1 — p over the 
whole range of p in the paramagnetic regime. The Hamil- 
tonian is H ~ ^ —TijSiSj where the sum is taken over 
all nearest neighbour bonds ij. 

For the discrete ± J interaction distribution the bond 
values are thus chosen according to 



P(J,,) = p6iJ,, - J) + (1 - p)SiJ,j + J) 



(1) 



We will set J = 1 and use when convenient inverse tem- 
peratures l3 — 1/T. The pure Ising model is recovered for 
p = I and the pure antiferromagnetic model for p = 0, 
while the standard strong disorder bimodal spin glass 
case corresponds to p ~ 1/2. 

In addition to the standard ferromagnetic and spin 
glass ordering transitions, on the RBIM phase diagram 
in the paramagnetic regime other physically significant 
lines can be operationally defined as functions of p. In 
addition to the Griffiths line (Tgrip) = T^) [l|,[3] and the 
Nishimori line Q there is the Fortuin-Kasteleyn (FK) 
transition line TpKip) 0-01 and the heat-bath damage 
spreading transition line Tds(p) [lol - fist . Early large scale 
relaxation measurements on the ISO in dimension three 



were interpreted |14l | in terms of a dynamic transition 
at the Griffiths temperature (i.e. Tq{p) = Tc), with a 



qualitative change in relaxation behavior of the autocor- 
relation function from non-exponential to exponential. 

It can alternatively be considered that Tp^ or T^j^ de- 
fines a dynamic transition. At TV if the FK cluster size di- 
verges so the standard cluster flipping algorithms 15 , 
break down. At T^s, the time scale for the "coupling 
from the past" equilibration criterion D{t) — )• diverges 



at large L [17 



18j so this criterion similarly breaks down. 
A priori it seems plausible that the two dynamic break- 
downs should occur at similar, if perhaps not identical, 
temperatures as has been conjectured by a number of 
authors 0,[ii,|23. 

In the following discussion it is shown that there are 
universal analytic relations for the intersection point be- 
tween Tpk{p) and the Nishimori line 01- With these 
exact values in hand together with the present accurate 
numerical damage spreading data a critical test can be 
made of the conjectured rule Tds{p) ^ Tpk{p)- This 
is shown to be a good approximation, almost exact in 
dimension three and accurate to a few percent in dimen- 
sions two, four and five. However the equality between 
the temperatures is not a general rule. 



TRANSITION DEFINITIONS 

Nishimori 0] has shown for the RBIM that, due to 
extra symmetries of the problem, a number of quantities 
may be computed exactly when the equality 



(1 -Pnl)/pnl = exp(-2^jvL) 



(2) 



holds; this condition defines the Nishimori line (NL), a 
line traversing the entire /3(p) phase diagram including 
both paramagnetic and ferromagnetic regimes. In par- 
ticular on the NL the internal energy of the system per 
bond (or edge) is 



Unl 



-[2pNL-l] = -tanh(/3. 



NL) 



(3) 



2 



This energy per bond has exactly the value that inde- 
pendent bonds would have at the same temperature. We 
will see below that indeed on the NL the positions of the 

and that bond 



satisfied bonds are uncorrelated 21 



We will refer to the temperature where the condition 



Pa{p,TFKr) ^ Pc 



(5) 



positions remain very close to random over a wide strong 
disorder regime around p = 1/2. In other words any- 
where on the NL the satisfied bonds are uncorrelated on 
average, so a fortiori the FK active bonds are distributed 
at random. 

The rule for the Fortuin-Kasteleyn (FK) transition line 
j^-Q) corresponding initially to the pure ferromagnets, is 
to first select for some particular equilibrium configura- 
tion the entire set of "satisfied" bonds ij where JijSiSj 
is positive. These bonds are then decimated at ran- 
dom leaving a fraction [1 — exp(— 2/T)] of "active" satis- 
fied bonds. The FK stochastic transition at Tpx occurs 
when the set of active bonds percolates through the lat- 
tice. From the way in which the bonds are laid down, 
in the general case this is correlated bond percolation, as 
opposed to the standard uncorrelated bond percolation 
with random bond occupation. Remarkably, for pure 
ferromagnets it can be proved that the FK transition 
Pfk{p = 1) coincides exactly with the Curie temperature 
Pc{p = 1) y. The Swendsen-Wang and Wolff cluster al- 
gorithms [la, [l^ speed up equilibration dramatically as 
long as the FK clusters remain of finite size. 

The same operational definition can also be used else- 
where in the phase diagram; in particular Tpk{p = 1/2) 
was estimated numerically for strong disorder RBIM on 
square, cubic, and triangular lattices It turned out 

that Tfk{p = 1/2) was in each case much higher than the 
standard spin glass ordering temperature Tg{p = 1/2). 
Even in the two-dimensional case where the glass tem- 
perature is zero, Tpk{p = 1/2) is similar to but lower 
than the Griffiths temperature Tgr = Tc{p = 1) [1, 0. 
The physical basis for the percolation transition temper- 
ature was explained in terms of "FK droplets" for the 
case of the fully frustrated lattice and general Potts q 
in Refs. 0, 24 1. Imaoka et al. [1^ estimated Tpk(j>) 



numerically over the entire range of p for the square and 
triangular RBIM lattices. 

Below we will introduce a further temperature TpKrip) 
closely related to Tpk{p)- By construction the fraction 
of satisfied bonds ris (p, T) in equilibrium is related to 
the equilibrium internal energy per bond U {p, T) through 
U{p,T) = 1 — 2ns{p, T) so the fraction of all bonds which 
are FK active bonds is just [l^ 

Pa{p,T) = (1 - C/(p,r))[l -exp(-2/T)]/2. (4) 

U(j),T) and hence Pa{p,T) can be readily measured nu- 
merically to high precision. It was noted [l^ that if the 
satisfied bond positions are assumed to be uncorrelated, 
the FK transition would occur at a temperature such that 
Pa{p,TpK) = Pc where Pc is the standard random bond 
percolation concentration for the lattice. 



is satisfied as TpKrip) {r standing for random). This 
conjecture led to estimates for TpKr{p = 1/2) on square 
and on cubic lattices in higher dimensions [Tot which 
were in good agreement with the numerical estimates for 
Tfk{p = 1/2) available at the time. 

It can be noticed that for the pure Ising square 
lattice ferromagnet at criticality U{f3c) = —1/2^^^ = 
-0.7071 . . . while - tanh(/3c) = -0.4121 . . . which is very 
different because the positions of the satisfied bonds are 
strongly correlated. However the critical FK concentra- 
tion of active bonds is Pa (p = 1, 7c) = 1/2 which is "acci- 
dentally" equal to Pc = 1/2 for this lattice. For the pure 
Ising ferromagnet on the triangular lattice at criticality 
on the other hand Pa(p = 1, Tc) 0.352208 . . . which is 
not quite equal to Pc — 0.347296 ... for this lattice. 

Finally, a heat-bath damage spreading transition, 
Tdsip) can be defined 13" 13 1- Heat-bath update rules 
are applied at fixed temperature to two initially non- 
identical spin configurations A{0) and P(0) (which are 
not necessarily equilibrated) of a given sample (mean- 
ing A and B have exactly the same sets of interactions 
Jij); the random number used in each subsequent single 
spin update step is the same for both configurations. At 
high enough temperatures, on annealing under this pro- 
cedure for sufficient time t, A{t) and B{t) will become 
identical; below the damage temperature T^s the "dam- 
age" {D{t) the Hamming distance between A{t) and B{t) 
divided by L'^) will stabilize for long times at a temper- 
ature dependent non-zero value. For the pure ferromag- 
net case Tdsip = 1) is equal to the Curie temperature 
Tc{p = 1) [13, [nj. For the strong disorder spin glass, 
just as Tfk{p = 1/2) is much higher than the glass tem- 
perature Tg{p = 1/2), Tds{p = 1/2) is also much higher 
than Tg{p = 1/2) |10l4l3| . The damage spreading tran- 
sition for a given coupling algorithm, in particular heat 
bath (HB), can be described as a "regular to chaotic" 
dynamic transition from the viewpoint of the "coupling 
from the past" approach. When the damage falls to zero 
after a sufficient anneal time it is a gua rantee that the 
system has been strictly equilibrated [3, [g^, [l^l • In sys- 
tems with frustration this guarantee breaks down when 
T < Tdsip) ^ so perfect equilibration in this sense for the 
regime near the critical temperature cannot be achieved. 

It should be noted that Tds depends on the updating 
protocol, with for instance Glauber updating giving very 
different results from HB updating. It turns out that 
even within the HB protocol the precise value obtained 
for Tds changes slightly depending on whether sequential 
or random updating is used. The results reported here 
are for random updating. With sequential updating the 
observed damage spreading temperature is of the order 
of 1% lower. 
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Early comparisons of numerical data from different 
groups indicated that Tds{p = 1/2) ~ Tpk{p = 1/2) 
for dimensions two and three 0, [l3 , and on this basis it 
has been conjectured by a number of authors 0, [3 US] 
that Tds(p) ~ Tpk{p) is a general rule defining a joint 
dynamic transition temperature above which relaxation 
is exponential in the long time limit and below which re- 
laxation is chaotic. It should be noted that in practice 
it is very hard to identify such a transition directly from 
autocorrelation function decay q{t) data. 

Recently Yamaguchi j23| focused attention on behavior 
on the NL. He provided conjectures suggesting that the 
intersection of the Tfk{p) hne and the NL would occur 
when 



the intersection. The present measurements show that 
in dimension two Tds{p) is in fact systematically lower 
than Tpk{p) by a few % over the entire range of p, ex- 
cept at and near the pure ferromagnet limit p = 1. At 
Pnl,fk for cubic models in dimensions three to five the 
conjecture holds to within about 0.1% in dimension three, 
to 2% in dimension four and to 4% in dimension five. 
The general relation Tds(p) ~ Tfk{p) is therefore simply 
a reasonably good approximation. In dimensions two, 
three and four the present data show that the damage 
spreading transition temperature like the FK transition 
temperature is very insensitive to p for a wide range of p 
around p = 1/2. 



PNL,FK = {l + Pc)/'2 



and 



NL.FK 



ln[(l + P,)/(l-P,)]/2 



(6) 



(7) 



where again Pc is the random percolation concentration 
for the lattice. It turns out that because of the analytic 
value for the energy Unl{p) on the NL, Eq. ([3]), the Ya- 
maguchi condition is strictly equivalent to the equality 
Tnl.fk ~ Tnl.fkt, meaning that this condition holds 
if the FK active bonds are distributed at random. The 
random bond conjecture for Pfk, Eq. (4) of Ref. is 
identical to Eq. (3.3)of Ref. 0. 

In fact Nishimori [21| years earlier had made the strict 
statement : "We have also proved independence of the 
local internal energy of different bonds, which indicates 
that the system effectively splits into uncorrelated sets 
of bonds on the [NL] in the phase diagram." In other 
words anywhere on the NL the satisfied bonds are un- 
correlated on average, so a fortiori the FK active bonds 
are distributed at random. Hence at the NL-FK inter- 
section concentration pnl,fk, we have Tfk = TpKrip), 
meaning that the conditions for this intersection conjec- 
tured by Yamaguchi, Eq. ^ and Eq. ([7]), hold exactly for 
any RBIM lattice if the FK transition is a random active 
bond percolation transition. Remarkably, a further con- 
sequence of Nishimori's statement is that this equality 
should hold for the mean over many samples indepen- 
dently of the sample size L. 

The numerical results below show that for dimension 
two, and presumably in other dimensions also, in ad- 
dition to this identity at pnl,fk, for a wide range of 
p around p ~ 1/2, Tfk{p) ^ TFKr{p) remains a very 
good approximation, as was conjectured for p = 1/2 in 
Ref. [19[. The active bonds are thus essentially uncorre- 
lated at Pfk{p) in this wide "high disorder" regime. It 
should be noted however that as a general rule the po- 
sitions of active bonds are correlated so Tfk{p) is not 
equal to Tfht^J))- 

Yamaguchi [20|, following Refs. 0, [3] for the case of 
p = 1/2, also conjectured that Tds{NL) = Tml.fk at 



FK NUMERICAL RESULTS IN DIMENSION 
TWO 

The two dimensional lattices studied were the square 
lattice and the triangular lattice. On both lattices the 
standard phase diagrams as functions of p are well es- 
tablished. The pure ferromagnetic Curie temperatures 
Tc{p = 1) arc exactly T^{sq) = -2/ln[2i/2 - 1] = 
2.26918... and T^itri) = 4/ln3 = 3.64096... respec- 
tively. When p is lowered from p = 1, the Curie tempera- 
ture T'c(p) drops gradually until a critical point is reached 
at pr on the Nishimori line where the Curie temperature 
tends suddenly to zero with weak re-entrant behavior 

[28-32|. Between the re-entrant regime and p = 1/2 
the system can be considered a "spin glass" but with no 
finite temperature ordering transition. The square lat- 
tice phase diagram for p < 1 /2 is the exact mirror image 
of the phase diagram for p > 1/2 with antiferromagnetic 
order taking the place of ferromagnetic order as the lat- 
tice is bipartite. For the triangular lattice on the other 
hand there is no finite temperature order below p = 1/2 
as the fully frustrated antiferromagnetic limit at p = is 
approached. 

The internal energy per bond U(j),T,L) at equilib- 
rium was estimated numerically. For small triangular 
lattices, 3 < L < 11, the energy as a function of temper- 
ature for each specific sample was calculated exactly (up 
to numerical precision) using transfer matrices. Once 
a sample is picked we split the lattice into L parts, 
each part corresponding to a transfer matrix Ai. Af- 
ter choosing a numerical value of T we evaluate the 
matrices and compute the trace of their product to ob- 
tain Z{p, T,L) = tr (Ai A2 • • • ,Al). This standard ap- 
proach is described in great detail for the Ising case in 
section 3 of [s^l and is of course easily adapted to the 
spin glass case. We thus computed (lnZ(p, T, i)) for 
p = 0, 0.02, 0.04, . . . , 0.50 and p = 0.51, . . . , 0.98, 0.99, 1 
and some 70-80 values of T. The number of samples 
ranged between 16384 for L = 3,4,5 and then down to 
only 1024 for L = 11. Taking the average dXnZ/dfi then 
provides us with {U{p,T, L)) . Values at {p,T) outside 
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the computed data grid were obtained through 3rd order 
interpolation. 

For larger L (for the square and triangular lattice) the 
Monte Carlo data were collected after equilibration us- 
ing standard Metropolis updating. At each {p,T), and 
for each sample, we collected at least a few million mea- 
surements of [/, and of course even more for moderate 
L. For the square lattices we used L = 16, 32, 64 on 64 
different samples with 0.5 < p < 1 in steps of 0.025. For 
the triangular lattice wc used L = 16, 32 on 32 differ- 
ent samples but with < p < 1 in steps of 0.025 (slightly 
denser at high and low p) . Some 50 values of T were used 
for both lattices. Again, intermediate values of (U) were 
obtained by interpolating the values at the (p, T) grid 
points. By interpolation of points derived from Eq. ^ 
and the consistency condition Eq. ([5]) the temperature 
TpKrip) where the fraction of FK active bonds Pa{p,T) 
is equal to Pc can be estimated to high precision. In 
Figures [1] and [2] the TpKr{p) values are compared to the 



directly measured Tpk{p) values [8|, |25| for the square 
and triangular lattices respectively. The damage spread- 
ing temperatures to be discussed later are also shown. 

The error in the Tpxr-estimates is found by first 
solving Equation ^ for each individual sample (using 
the interpolation we mentioned above) which gives us 
a sample-to-sample standard deviation from which we 
obtain the standard error of the average-sample TpKr- 
Once the standard error is established we are free to plot 
a smoothed version (from fitting a polynomial of high 
degree) of TpKr versus p. Needless to say we have not 
estimated the error in [/, and hence Tpxr, for each in- 
dividual sample, but we simply assume this error should 
be reflected in the sample-to-sample variation. 

In dimension two, Pc = 1/2 for the square lattice 
and Pc = 2sin(7r/18) = 0.347296355... for the tri- 
angular lattice. The exact NL-FK intersection values 
from Eqs. dH) and ([7]) are then : pnl,fk = 3/4, 
pNL,FK = ln(3)/2 = 0.5493063 ... for the square lattice, 
and pnl,fk ~ 0.6736 . . ., (Unl.fk — 0.362371 ... for the 
triangular lattice. The estimated Tpk{p) and Tp^rip) 
lines run directly through these exact NL-FK intersec- 
tion points as they should. Figures [1] and [2l 

The measured equilibrium energy per bond is slightly 
higher than the random bond energy — tanh(/3) for /? > 
f3pKr and slightly lower for (3 < jSpKr- 

To within the high numerical precision the estimated 
PpKriPNL.FK, L) at the NL-FK intersection point is in- 
dependent of L down to i = 3, Figures [3] and ID The ob- 
served absence of finite scaling corrections at the NL also 
follows from Nishimori's general statement quoted above 
PH . For p- values to the left and right of pnl.fk there 
are weak finite size effects of opposite signs. This absence 
of mean finite size scaling deviations arises because TpKr 
depends only on the energy. Other parameters can still 
show finite size scaling deviations |34l | . The large size lim- 
its for Tpk{p) and TpKrip) remain equal to each other 
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FIG. 1: (Color online) The square lattice transitions, tem- 
peratures T as functions of ferromagnetic interaction concen- 
trations p. The pink dashed line on the right indicates the 
ferromagnetic Curie temperatures, data taken from Ref. [s^ . 
The green solid traversal line is the exact Nishimori line (NL). 
Red circles : FK transitions Tfk, Ref. [1^. Blue squares : 
random active bond percolation line Tfkt- Black diamonds 
: heat bath damage spreading transition Tds- The errors on 
Tfk and Tds are of the size of the points. Large black square 
: exact intersection point Tnl,fk- 



to within the numerical precision of the Tpxip) points 
[2^ over the strong disorder range of p extending from 
p = 1/2 to pnl,fk in the square lattice, and from p ~ 0.3 
to pnl,fk in the triangular lattice. Over these ranges of 
p the satisfied bonds at TpK are very close to being un- 
correlated; for instance at p = 1/2 on the square lattice 
—U{p,/3pKr)/ ta.nh{f3pKr) ~ 0.9709 which remains near 
to the uncorrelated value of 1. The further randomness 
introduced by the FK decimation is sufficient to render 
the active bond positions essentially uncorrelated. The 
Tpk{p) [23(or TpKr{p) which has been measured here 
with higher precision) is rather insensitive to p within 
this range, so that the measured Tfk(1/2) is similar to 
the exact Tml,fk on both lattices. For both lattices in 
the more strongly ferromagnetic ranges p > pnl,fk, and 
in the range p < 0.3 near to the fully frustrated p = 
limit in the triangular lattice, Tpk{p) > Tpxrip)- At 
these concentrations the satisfied bonds are significantly 
correlated at Tpk{p) and the FK decimation does not 
sufficiently compensate the correlations so as to produce 
randomness among the active bonds. 
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FIG. 2: (Color online) The triangular lattice transitions, tem- 
peratures T as functions of ferromagnetic interaction concen- 
trations p. The pink dashed line on the right indicates the 
ferromagnetic Curie temperatures. The green solid traversal 
line is the exact Nishimori line (NL). The red dashed vertical 
line in the centre incates p = 0.5. Red circles : FK tran- 
sitions Tfx, Ref. [2^. Blue squares ; random active bond 
percolation line Tfkt- Black diamonds : heat bath damage 
spreading transition T^a- The errors on TpK and T^a are of 
the size of the points. Large black square : exact intersection 
point Tnl,fk- 




FIG. 3: (Color online) The size effect on triangular lattice 
temperatures Tfkt at which the fraction of FK sites is equal 
to the random site percolation concentration for the lattice 
Pc, Eq. ((5)1. Temperatures Tfkt as functions of ferromagnetic 
interaction concentrations p for lattice sizes L from 3 to 11. 
L = 3,4,5,11 black, red, green, blue; from top to bottom 
on the right hand side and from bottom to top on the left 
hand side. The black square is the exact Tml.fk intersection 
Eqs. Q and ©. 



FIG. 4: (Color online) As for Fig. [3l a closeup of the re- 
gion around the Tnl.fk intersection (square). L = 3, 4, 5, 11 
(black square, red circle, green triangle, blue diamond; from 
top to bottom on the right hand side and from bottom to top 
on the left hand side). The error bars are the standard error 
when measuring Tfkt after solving Eqs. (jj) and © and cor- 
respond to 16384 samples for L = 3, 4, 5 and 1024 for L — 11. 
The curves are 10th degree polynomials fitted to the whole 
range from p = to p = 1. 



DAMAGE SPREADING NUMERICAL RESULTS 
IN DIMENSION TWO 



For dimension two the critical temperatures T^s (p) be- 
low which the long time large L heat bath damage spread- 
ing D(j), T, t) tends to a non-zero value are also shown in 
Figures [1] and [H The values given are obtained from 
extrapolating data at increasing sizes, Figures [5] and [Bl 
to infinite size. On the basis of data for square and cu- 
bic lattices at p = 1/2 it was conjectured earlier that 
TpK ip) ~ Tds (p) 0, Hi E^l as is the case in the pure fcr- 
romagnet p = 1 limit. However the present data demon- 
strate that while the curves for Tds{p) and Tpk{p) he 
near together and are of very similar shape, in particular 
both being almost independent of p for the range near 
p = 1/2, it is clear that Tds{p) < Tfk{p) except when p 
tends to the ferromagnetic limit p = 1. So unfortunately 
it is not possible to define a unique joint "dynamic tran- 
sition temperature" . We will see that this is true also in 
higher dimensions. 

The present accurate critical value for the strong dis- 
order square lattice T(is(l/2) = 1.69(2) is in good agree- 
ment with the earlier value Tdsi)^/'^) = 1-70 [l3| and is 
very close to the " regular to chaotic" dynamic transition 
temperature estimated for the same lattice in Ref. [3| 
from the divergence of the coupling time with LF' , where 
the data indicate a transition at Tds{^/'2) ^ 1.72. 
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FIG. 5: (Color online) The damage spreading D{l3, L) for 
the triangular lattice ai p — 0.6736. The lattice sizes are 
L = 48, 64, 96, 128, 256, 384, 512, 1024 (black square, red cir- 
cle, green triangle, blue inverted triangle, cyan diamond, yel- 
low left triangle, pink right triangle, brown star. Right to 
left in the region where the curves separate). Shown error 
bars are the standard error. The exact FK transition inverse 
temperature is /3 = 0.362371 ... so clearly distinct from the 
limiting damage spreading inverse temperature. 



DIMENSIONS THREE, FOUR AND FIVE 

It is numerically much more demanding to esti- 
mate Tfk (p) precisely through direct measurements 
(as in [251 ), in particular allowing for finite size ef- 
fects, than to estimate Tds{p) to the same level of ac- 
curacy. For cubic lattices in dimensions three, four 
and five the random bond critical concentrations Pc 
though not exact have been estimated to very high pre- 
cision, Pc(3) = 0.2488126(5), Pc(4) = 0.1601310(10) 
and Pc(5) = 0.11811718(3) respectively [H, E^. Us- 
ing the exact NL-FK intersection point expressions, 
Eqs. and ([7]), one thus has pnl,fk = 0.6244 and 
Pnl.fk = 0.2541466 in dimension three, pnl,fk = 
0.58006 and Pnl.fk = 0.161521 in dimension four, and 
PNL,FK = 0.55906 and I3nl,fk = 0.118671 in dimen- 
sion five. Hence as pnl,fk and I3fk{pnl,fk) can be 
taken as known almost exactly, it is sufficient to esti- 
mate l3ds{PNL,FK) numerically (allowing carefully for fi- 
nite size corrections) to obtain an accurate estimate of 
the ratio Tds(jpNL,FK) /Tfk{.Pnl.fk) for each dimension. 

The equilibrium damage D{L,I3) was measured for 
given L as a function of /3, and the results were extrap- 
olated to obtain an estimate of the /3ds value at which 
D{oo,(3) falls to zero. Figures [5] and [6l There is a clear 
envelope curve for all L in each case, with the data points 
leaving the curve later and later as L increases. There 
are various ways to extrapolate to infinite L. One cffi- 



FIG. 6: (Color online) The cubic lattice equilibrium damage 
spreading D{(3, L) as a function of size and inverse tempera- 
ture for p = 0.6244. Sizes : i = 12, 16, 24, 32, 48, 64 (black 
square, red circle, blue triangle, green inverted triangle, cyan 
diamond, pink star ; right to left in the region where the 
curves separate). Shown error bars are the standard error. 
The extrapolated intersection with the /3 axis gives the in- 
finite size critical Pds- In this case fids = 0.25432(15) can 
hardly be distinguished from the exact FK transition inverse 
temperature /3fa' ~ 0.25414 . . ., red arrow. 

cient method is to fit the envelope points assuming that 
near /3ds the behavior follows D{/3) = A{/3 — fids)^ and 
adjusting B to obtain a straight line. We have also used 
a simple scaling formula PdsiL) = Pds + C L^^ to ver- 
ify. Of course, a certain statistical uncertainty will enter 
depending on which L to include in the fitting process 
which we take into account in the final error estimate. 

The infinite L damage spreading temperatures 
Tds{pNL,FK) were estimated for the central values of 
Pnl,fk in the three dimensions : Tds{pNL,FK) = 
3.932(2) in dimension three, Tds[pNL,FK) = 6.057(10) in 
dimension four, and Tds{pNL,FK) = 8.13(1) in dimension 
five. The observed ratios TFif(pArL,Fif)/7ds(pArL, fx) are 
equal to 1.0005(5), 1.022(2) and 1.036(10) respectively. 
In dimension three Tds{pNL,FK) is indistinguishable from 
Tfk{pnl,fk), while in dimensions four and five the val- 
ues appear tantalizingly close but not identical. We have 
no explanation for the striking similarity between the 
two temperatures for the particular case of dimension 
three. Adding a little more detail for this particular case, 
the individual j3ds{L) were estimated to an accuracy of 
±0.0001 for L = 48, 64, ±0.00025 for L = 32, ±0.0005 
for L = 16, 20, 24 and ±0.001 for L = 12. Using the ex- 
ponent A = 2.05 gives a projected (3ds that depends only 
to a very small degree on which L are included in the fit 
(though we always include L = 48, 64). We receive a me- 
dian value of /3ds — 0.25432 and a standard deviation of 
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0.00015 giving us Tds = 3.932(2). Similar methods were 
used to estimate Tds for the other dimensions. 

As in dimension two, the curve for T^sip) is ahuost flat 
over a wide range of p around p = 1/2 for dimensions 
three and four, Figures [7] and H] (the case of dimension 
five was not studied). These data are for fixed L (L = 32 
and L = 16 respectively) so the Tds{p) values are not 
quite equivalent to the values quoted for Tds at Pnl,fk- 

Judging from the behavior observed in dimension 
two, it seems very plausible to assume that the ratio 
TpK (p) /Tds (p) is always practically independent of p in 
the strong disorder regime 1/2 < p < pnl,fk- 

It can be noted that numerical data for the time depen- 
dence of the autocorrelation function g(t) at all temper- 
atures in two dimensional fully frustrated systems have 
been interpreted in terms of exponential relaxation with 
logarithmic corrections due to vortex-vortex interactions 
[37|. The results indicate that in these systems (which 
have well established FK and damage transitions) there 
is no dynamic critical temperature in the Ogielski sense 

0- 
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FIG. 7: The damage spreading critical temperature Tds{p) 
for a 3d cubic lattice size L = 32. The error bars correspond 
mainly to the residual uncertainty in the extrapolation to the 
intersection with the /? axis. 



DAMAGE CLUSTERS 

It is of interest to examine in some detail the damage 
clusters in dimension two where they can be readily vi- 
sualized. Keeping (3 > Pds and letting t ~¥ oo (or at least 
allowing for equilibration) how does the damage actually 
spread? Note that we define a cluster as a connected 
component in the lattice induced by the damaged sites, 
i.e. a maximal set of damaged sites such that there is a 
lattice path (made up of horizontal and vertical steps) be- 
tween each pair. The sum of the individual cluster sizes 
is thus the damage without normalization, i.e. L'^ D{f3). 
We have collected data on cluster sizes for d ~ 2, 3, 4, 5. 
For d = 2 we have only used the square lattice, not the 
triangular lattice. 

What we see in a square lattice snapshot at (3 near (3ds 
and at any arbitrary fixed time is not, as in the case for 
percolation (including FK percolation) or random graphs 
near Pc, the formation of a giant cluster, but rather a 
large number of very small clusters. At subsequent times 
the clusters evolve and flutter through the lattice. An 
example is shown in Fig. [3] 

The measurements take place as above for the damage 
spreading; starting from random (infinite temperature) 
spin configurations, at each time step we update L'^ ran- 
domly selected sites, then search for all clusters and their 
sizes. On a measurement we collect the number of clus- 
ters, the average cluster size, the size of the largest cluster 
and the size of a randomly selected cluster. As before, 
we collect 10^ measurements from each sample after dis- 
carding between 50000 and 250000 time steps (depending 
on L) to allow for equilibration of the damage spreading. 
We then average over the time steps to get the average 
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FIG. 8: The damage spreading critical temperature Tds{p) 
for a 4d cubic lattice size L = 16. The error bars correspond 
mainly to the residual uncertainty in the extrapolation to the 
intersection with the /? axis. 

for a particular sample. The data are then averaged over 
the samples; we have only used eight samples in each case 
for the cluster measurements (for the damage speading 
we used between 8 and 128 samples depending on lattice 
size). In all cases we have set p = pnl,fk- 

First we discuss the expected number of clusters, which 
we denote nc(/?, L). In Figure fTOl we plot nc(/3, L)/L^ ver- 
sus P for the simple cubic lattice. Equivalent behavior 
is found also for d = 2,4,5. Note that for the square 
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FIG. 9: An illustrative snapshot of an instantaneous damage 
site configuration for an L = 128 square lattice at p = 0.75 
and temperature T — 1.61, just below the damage spreading 
transition. 



lattice all lattice sizes agree on a global maximum prob- 
ability located at /3max = 0.656(1), where nc(/3max,i) ^ 
0.0130(1)2.2. For d = 3 we obtain ^^ax = 0.313(1) and 
nc(/3max,i) ~ 0.0227(1) for d = 4 we get /3„,ax = 
0.211(1) and nc(/3max,i) ~ 0.0200(1) L", while d = 5 
gave Anax = 0.177(1) and ne(Anax,i) - 0.0308(1) L5. 
Thus the damage is always distributed over 0{L'^) clus- 
ters. 

With so many clusters the average cluster at each 
time step must be rather small. We measure the num- 
ber of damaged sites D L'^ and the number of clusters 
ric, giving us the mean cluster size at each time step. 
The normalized time average (and then the sample aver- 
age, though the sample variation is very small), is then 
s(/3,L) = {DL'^/tIc)- In Figure[Tl]we plot s{/3,L) versus 
/3, again for d ~ 3. The curve clearly suggests a posi- 
tive right limit at /3ds- To estimate this limit we fit a 
simple expression + ai exp(a2 /3) to the points result- 
ing in s{l3ds) = lim^^^j- limi^oo s(/3, L) = 3.46(1). We 
estimated the right limits corresponding to d = 2,3,4,5 
to be respectively 9.00(2), 3.46(1), 2.76(1) and 1.95(1) 
sites at /3ds- One could alternatively define s as the av- 
erage damage divided by the average number of clus- 
ters. This is not strictly the same as our present def- 
inition (mean ratio versus ratio of the means) but the 
difference is of course vanishingly small here. For exam- 
ple, the maximum in Figure [10] at /^max = 0.313 gives 
Tie ~ 0.0227 and the number of damaged sites, see 
Figureini is 0.139(1) L^. Hence the average cluster size is 
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FIG. 10: (Color online) The number of clusters nc{l3,L) on 
the 3d cubic lattice normalized by at p = pnl,fk- Sizes 
L = 12, 16, 24, 32, 48 shown as black squares, red circles, blue 
triangles, green inverted triangles, pink left triangles. The 
errors are smaller than the size of the points. 
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FIG. 11: (Color online) The mean cluster size s{f3,L) on the 
3d cubic lattice at p = pnl,fk- Sizes L — 12,16,24,32,48 
shown as black squares, red circles, blue triangles, green in- 
verted triangles, pink left triangles. The errors are smaller 
than the size of the points. 



roughly 6.1 which matches Figure [TT] 

We end this section with a final remark concerning the 
distribution of cluster sizes. At each time step we here 
pick one cluster uniformly at random and measure its 
size. In Figure [T^ we show a set of size distributions 
(or density functions) as log Pr vs log s for a range of f3 
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FIG. 12: (Color online) The distribution of cluster sizes on 
the 3d cubic lattice for L = 12 at inverse temperatures 
P — 0.295, 0.335, 0.370, 0.500 (black squares, red circles, green 
triangles, and blue inverted triangles). The errors are of the 
size of the points except for very small p where they can be 
judged by the scatter (note the logarithmic scales). 



for d = 3, L = 12. They all have a high value at size 1 
(many isolated sites) and then drop quickly. They can 
be expected to behave like this since the mean size is 
between 4 and 70 in this temperature range. However, 
for f3 > /3max (or thereabout) the distributions show a 
second, and rather wide, maximum located in the neigh- 
borhood of the total damage D{/3, L) . This can be under- 
stood as follows. Suppose for the sake of argument that 
the damaged sites are distributed at random in space. 
(This is only approximate as correlations between dam- 
age sites should be allowed for). Then when D(/3) exceeds 
the site percolation concentration there will exist a sin- 
gle "giant" percolating cluster of damaged sites together 
with residual small clusters. As D{(3) increases further 
the percolating cluster will contain almost all the dam- 
aged sites. On this criterion the peak should appear in 
the distribution at D{I3) - 0.59,0.31,0.20 and 0.14 in 
dimensions 2, 3, 4 and 5 respectively which gives an indi- 
cation in rough agreement with the data. Indeed in the 
case of the square the peak never appears which is consis- 
tent with the fact that D{j3) never approaches 0.59. For 
the dimensions where the peak does appear the global 
parameters, in particular D{p), increase smoothly with 
(3 and show no sign of any critical behavior as the giant 
cluster forms. 



CONCLUSION 

The exact values of the coordinates of the intersec- 
tion point where the Fortuin-Kasteleyn transition line 
crosses the Nishimori line can be derived for an RBIM 
from the analytic condition that satisfied bonds arc un- 
corrclatcd on the NL |2lJ . On any lattice this leads to the 
exact expressions pnt,.fk = (1 + Pc)l^ and Tj^l^pk ~ 
2/ln[(l -(-Pc)/(1 --Pc)j EliQi where Pc is the standard 
random bond percolation concentration for the particular 
lattice. For lattices in dimension two (and probably in 
higher dimensions also) the uncorrelated bond condition 
remains a very good approximation at the FK transition 
temperature over a wide strong-disorder region spanning 
p=l/2. 

In pure ferromagncts Tpk{p = 1) = Tds{p = 1) = 
Tc{p = 1). The conjectured equivalence for the RBIM 
between the FK transition temperature and the heat 
bath damage spreading temperature , Tpk{p) ^ Tds{p) 
0, ESl : separating an exponential from a chaotic dy- 
namic regime has been tested at p = pnl.fk on simple 
cubic lattices in dimensions two, three, four and five. It 
holds to within 0.1% in dimension three, to within 2% in 
dimension four, and to 3% in dimension five. The equiv- 
alence appears always to be a good approximation. We 
have no explanation to propose for the quasi-cquality in 
the case of dimension three. 

For the square and triangle lattices the difference is 
larger : 5% and 16% respectively. 

The FK transition in the strong-disorder regime close 
to p ~ 1/2 can thus be said to be well understood. How- 
ever the basic physical condition determining the the 
damage spreading transition temperature, which plays an 
important role in limiting perfect equilibration in RBIMs, 
and its proximity to the FK transition, remains unclear. 
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